MULTIFOCAL OPTICAL DEVICE DESIGN 



Technical Field 

The present invention relates to computer-implemented methods and systems for 
designing optical devices, and particularly concerns the design of multifocal devices such 
as progressive ophthalmic lenses. 

Background 

Multifocal optical devices find utility in a number of technologies. For example, 
bifocal and trifocal lenses treat presbyopia caused by loss of accommodative power of the 
eye. Progressive addition lenses for this purpose include a far-view zone with a low 
power for distance vision and a near-view zone with higher power for reading and other 
close work. The power increases progressively and smoothly in an intermediate zone or 
corridor between the far and near zones. Normally, the near zone lies near the bottom 
center of the lens, below the far zone; however, other arrangements satisfy specialized 
purposes. 

Lens design generates a specification of surface points that a grinder can produce 
fi-om a lens blank. Modem lens grinders are numerically controlled machines that abrade 
the front and/or back surfaces of the blank according to design data specifying various 
points on the surface. Some lenses may alternatively be produced by plastic casting or 
other techniques. Both the firont and back surfaces of a lens may have magnifying power; 
the effective power is approximately the difference between the powers of the two 
surfaces. Power at any point is related to the sum of a surface's curvatures in two 
perpendicular directions. 

Astigmatism is normally an undesirable property of lenses that causes distortion 
and color fringing. Astigmatism, sometimes called "cylinder," is related to the difference 
between the lens curvatures in different directions. While a single power lens can 
theoretically have zero astigmatism everywhere, a lens having different powers at 



different locations necessarily has at least some astigmatism. A major goal in designing 
progressive lenses is to reduce astigmatism, especially in the region of the far, near, and 
intermediate zones. Secondary zones outside this region are less critical, although 
astigmatism reduction there is still important. 
5 Computer-aided design of progressive lens and similar devices have employed 

both direct and indirect methods. A direct method first prescribes a selected distribution 
of curvatures on a vertical meridian line that follows the eye's up-down motion. 
Prescribed horizontal curves across the meridian are chosen to have the desired 
curvatures where they cross the meridian. Direct methods usually result in a high 
10 astigmatism level with little control over the astigmatism in different areas of the lens 

surface. 

Indirect methods balance a desired power distribution with the unavoidable 
astigmatism produced. A variational approach iterates toward the solution of a full non- 
linear fourth-order Euler-Lagrange partial differential equation. Another approach 

15 attempts to minimize a defined cost function with an optimization algorithm. The cost 

function attains a low value when the desired lens properties are approached. One such 
method specifies a desired power distribution over an entire lens surface, then constructs 
a surface close to this distribution while optimizing to achieve low weighted amounts of 
astigmatism. Free-energy analogies and variational approaches have also been 

20 investigated. Indirect methods involve the solution of fourth-order nonlinear elliptic 

equations, which involve very large amounts of computation — ^up to several hours for a 
single lens design. 

A lens design may deliberately include an astigmatic or cylinder component to 
correct for astigmatism in the wearer's eyes. Integrating the magnifying-power correction 
25 and the astigmatism or cylinder correction in the same surface of a lens or other optical 

device is desirable. Single-power lens conventionally fabricate both corrections in the 
same surface of a toric (or atoric) lens. However, progressive design techniques have 
heretofore not integrated astigmatism correction, and fabricating the two corrections 
sequentially on the same surface would unacceptably thin the lens or device. 
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Conventional methods therefore fabricate the progressive component on one surface 
(usually the front), and the astigmatism correction on the other sxirface. Although this 
makes the two corrections approximately independent of each other, some offset errors 
still exist, and structural weakness may still occur, especially in currently fashionable 
5 thin, light ophthalmic lenses. 

Summary of the Invention 

We offer an indirect way to design progressive optical devices that reduces 
computation time for accurate powers and small undesired astigmatism. Input parameters 
1 0 comprise a prescribed power distribution, a weight function for power deviation, and one 

for astigmatism. We solve a linearized variational equation as perturbations over a base 
surface, and may express the solution in terms of superposed low-degree splines. 

Another aspect of the invention permits progressive-power and astigmatism 
corrections in the same optical surface by incorporating a prescribed cylinder into the 
15 base surface. 

Drawing 

Figs. 1 A-E depict an example of a progressive lens design according to an 
embodiment of the invention. 
20 Fig. 2 is a flowchart for an example of a method for designing lenses such as that 

of Fig. 1. 

Fig. 3 shows the partition of a design field into zones. 
Fig. 4 illustrates a grid for finite element analysis. 
Figs. 5A-B show a spline function. 
25 Figs. 6A-B depict an example of a progressive lens design incorporating 

astigmatism correction. 

Fig. 7 is a block diagram of a system for hosting the invention. 

Detailed Description 
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Progressive Correction 
Fig 1 shows an example ophthalmic progressive addition lens (PAL) 100 inside an 
edge 1 10. This lens typifies designs made with our procedure, although other optical 
devices can be produced in a similar manner. The coordinate axes denote millimeters 
5 from the center of the lens. Fig. 1 A includes contour lines 120 having constant optical 

power. The mean power of an optical surface at a point x,y is proportional to the sum of 
the curvatures in two perpendicular directions at that point; the mean curvature is defined 
as H{x,y) = {k^(x,y)'^k2(x,y))/2. In the example lens 100, the add power is 0.00 diopter in 
the far-view upper part of the lens. The near-view add power is 2.0 diopters in the 

10 bottom center of the lens. ("Add power" is the power that must be added to the far-view 

power to obtain the desired near- view power; in this example, the base or far- view power 
is 4.94.) A vertical corridor between the near-view and far-view portions has powers 
intermediate the two specified powers. 

Fig. IB includes contours 130 of constant astigmatism, sometimes called cylinder. 

1 5 Astigmatism is proportional to the curvature difference between the two directions at the 

point , A(x,y) = \ki{x,y)~k2{x,y)l Unless astigmatism is deliberately designed in to correct 
an astigmatic condition, it should be lowest in the far- view, near- view , and corridor 
portions. Astigmatism in other portions of the lens is not as critical, although a good 
design attempts to reduce it throughout the lens surface. 

20 Fig. 2 is a flowchart showing an embodiment of a method 200 for designing a 

multifocal optical device. In this embodiment, the device is a progressive ophthalmic 
lens. 

Block 210 accesses a predefined design field. The field may be partitioned into 
zones for the purpose of assigning weights indicating the importance of various locations 
25 in the field. Fig. 3 shows a representative field 300 over which method 200 produces 

output data. Most lenses have a generally circular shape, and the field can be defined in a 
circular or other overall shape. However, we have found that a square field allows us to 
simplify solution of the equations. Dashed line 310 indicates the edge of a typical lens 
within the field. Output data outside this edge is simply discarded. The axes of field 300 
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indicate coordinates in millimeters away from the lens center. A primary region 320 
defines areas of higher importance, in which the astigmatism should be low. Region 320 
includes a zone 321 roughly corresponding generally to the base-power far-view area of 
an ophthalmic lens, a zone 322 in the area of the add-power near-view zone, and an 

5 intermediate zone 323 forming a corridor of progressively increasing power between 

zones 321 and 322. A secondary region 330 defines field areas of lesser importance, such 
as outlying zones 33 1-334. Region 330 at least partially surrounds region 320. More or 
fewer zones and differently shaped zones are also useful. 

Block 220 receives a data from the designer to specify parameters for a particular 

10 design. Some of these concern properties of the lens. For example, the designer specifies 

the overall size of the lens and the index of refraction n of the lens material. Typical 
diameters are about 60 to 80 mm. For common optical glasses used in PALs, l.S^n^lJ 
is a typical approximate range. An add power for near vision is specified in diopters, 
commonly 0.75 to 3.5 in increments of 0.25. Some parameters concern characteristics of 

1 5 the design, such as a base radius of the lens. For ophthalmic applications, this radius 

usually lies in the approximate range of 65-350mm, smaller radii being usual for high 
prescriptions (e.g., base power +3), and larger radii for lower prescriptions (e.g., - 10). A 
designer commonly selects appropriate sizes, front-surface powers and back-surface 
powers from standard charts for desired prescription ranges. The grid size of field 300 

20 may be allowed to vary as an input parameter as explained below. Although there is no 

inherent limitation on grid spacing, smaller spacings have higher resolution, but larger 
ones favor computational efficiency; acceptable compromises typically make each grid 
square or patch approximately l-4mm. 

The designer selects a desired or prescribed power distribution function P{x,y) 

25 over design field 300. The value of this function approaches the far-vision power in zone 

321 and the near- vision power in zone 322. The designer may encode into the power 
function such specifications as the length and angle of the corridor between the far-vision 
and near-vision zones, and the point at which the power begins to change. This function 
may be specified in a relatively arbitrary manner over the lens surface. It is often 
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convenient to specify it as a piecewise constant function having different values in each 
zone. Because of the sharp discontinuities of such a function, it may be smoothed, for 
example by convolution with a Gaussian function or a spline, if desired. Fig, IC graphs a 
power-distribution function 140 for example lens 100. 
5 The designer also assigns two weight functions over the area of field 300. A 

designer selects their magnitude distributions in the x and y directions of field 300 for a 
particular design. The weight functions may be in any form. It is often convenient to 
assign them as piecewise-constant functions, with or without smoothing, conforming to 
the boundaries of the zones defined above. An astigmatism weight function a{x,y) is 

10 large where the designer wishes small astigmatism. A power-deviation weight function 

I3{x,y) is large where the designer wishes low deviation from the power prescribed by 
Pix,y), In this example, both weight functions are positive-valued, in order to guarantee 
existence and uniqueness of the design solution. A designer typically wishes to overlap 
areas of low astigmatism with areas having low power deviation. Providing two 

15 independent weight functions allows freedom to assign these areas in any desired manner 

with respect to each other for different applications and preferences, and allows choices 
as to power accuracy and astigmatism to be made separately. A lens configuration 
depends strongly upon these two functions. The high computation speed of method 200 
allows a user to pursue many variations of a design. Figs. ID and IE graph sample a and 

20 )3 weight functions 141 and 142 for lens 100. 

Blocks 230 comprehend operations for generating output data representing the 
surface of an optical device such as progressive or multifocal lens 100. For a surface 
defined parametrically as u{x,y), we have selected a measure of the siuface quality 
defined as: 

25 /(«)= J |a(x,>')(^)\;ff(x,;.)(^-P(;c,>.))'|^rfy 

where the integral is taken over the entire design field 300 and P is the prescribed power at 
various points in the design field. (All integrals of unspecified range herein are understood to 
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be taken over the entire design field 300.) An overall lens design depends significantly upon 
the designer's choices of the weight functions amd /), This gives the designer fii-eedom to 
select among many options for different applications and preferences, and allows choices as to 
power accuracy and astigmatism to be made independently of each other. 
5 In this context, optimizing the design means minimizing the value of the surface- 

quality measure, which may also be thought of as a cost function. Its global minimum 
specifies the height of a surface z=u(x,y) within the design field 310 that has nearly correct 
powers and low astigmatism in the areas selected as important by the weight functions. For 
the present purpose, we have modified the formulation of this measure so as to express it in 
10 terms of the mean curvature H(x,y)y as defined above, and the Gaussian curvature K{x,y) = 

ki{x,y)>^k2{x,y). The quality measure then becomes: 

\la(H{x,yf - Kix.y))^ Wix.y)- P{x,y)f]dxdy 

The explicit forms of H and K, in terms of uix^y), then become: 

(1 + u])u^ - lu.UyU^ + (1 + ul)u^ 



15 where the subscript notation signifies partial differentiation with respect to the subscripted 

variable. From the variational calculus, the necessary condition for a minimum in I{u) is: 



dt 



20 



for all infinitely differentiable functions S{x,y) having compact support. This condition 
translates to the following form: 

J{ 2((a + m-miH.A + + H„J„ + H^J^ + H,J^) 
- a(K.S, + K^^S^ + KJ^ + K,J^ + KJ^) }dxdy = 0 
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If a, and P are at least twice differentiable, then integration by parts twice yields the fourth- 
order partial differential equation: 

l{{aH + PH - pP)H^^ \ - l{{aH + - pP)H^^ \ 
+ 2{{aH + - PP)H^^ )^ + 2((a// + pH - )^ 
+ 2((ai/ + ;8//-^i>)//„J^ = 0 



This is a form of the Euler-Lagrange equation corresponding to this variational problem. This 
5 Euler-Lagrange form is nonlinear, requiring excessive amounts of computational power and 

time for solution. A finite-element solution of acceptable spatial resolution might require 
several hours on the types of computers now economically feasible for use in this type of 
application. However, we have constructed a simplification that greatly reduces the 
computational load while retaining accuracy and resolution. The surfaces for most 

10 ophthalmic prescriptions, and for many other kinds of practical lenses, are almost spherical. 

We therefore treat the surface u{x,y) as a perturbation v{x,y) from a specified spherical surface 
u^{x,y)-^{B^'0^'y^) of radius jR, and we treat the gradients and of the perturbation as small. 
(Although the equation could alternatively be linearized about a plane or other fixed surface, 
this embodiment employs spherical linearization as preferable for many ophthalmic 

15 applications.) These modifications have the effect of linearizing the Euler-Lagrange equation, 

which allows us to express H and K as: 

1 (1 + ^lx>yy ' 2%^o,V^ + (1 + Uly)v^ 
^uW = 7 + 9/1 . „2 . 2 x3/2 + ^1 



where and E2 are small error terms. We select a domain where x^+y <i?^ where R is the 
radius of the spherical surface Mq. This radius may be taken as the base radius specified above 
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as a parameter. Then, if (v^vj«l and if\v^v^vj is bounded from above, E, and become 
insignificant. 

We then recast the variational problem with these simplified forms, using a constant 
spherical surface of radius R that is selected by the designer or chosen in some other way. 
The above necessary condition for a minimum of I(u) then becomes: 



f [2(a + m^MH^iS)- aK^iS,v)]dxdy 



The terms and K„g, which no longer represent true curvatures, are defined as: 



(1 + ul)v^ - 2u,,Uo^v^ + (1 + uly)v^ 
2(1+Mo'. + 0' 



Tj / \ V- "vxjyy vx 



S V +S V -2d V 

XX yy "yy xx xy xy 



V (X .A " " ^ 



10 



Here again, S(x,y) is an arbitrary test function expressing the variation. This formulation 
demonstrates that is linear, and K^q is bilinear and symmetric. We then defined two 
fimctions B{d,v) and L{6) as: 

B{5,v)= J [2ia+fi)H^(S)H^iv)-aK^iS,v)]dxdy 



aK(d,u,)+2\pP- 



R 



dxdy 



where 5 is a bilinear form depending on 6 and the perturbation function v, and where L is a 
1 5 linear functional of 6 alone. We then pose the minimum cost as a variational form 

B{6,v)-L{S), for all 6 in the relevant subspace of functions having the desired properties. 

Block 231 of method 200 solves the linearized partial differential equation 
B{S,v)-L(S) as a boundary- value problem for the parameters input in block 220. The solution 
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v(x,y) represents the perturbation or difference in height from the specified spherical surface 
Mo(x,3/) over the entire design field 300. 

A number of general finite elements are known for solving the present type of fourth- 
order equations on a digital computer, such as the Argyris element and the Hsieh-Clough- 
Tocher (HCT) element. However, these are computationally intensive. The Argyris triangle, 
for example, employs polynomials of degree five. Although the HCT element uses only cubic 
polynomials, each triangle must be divided up into three smaller triangles, each computed 
individually. As noted above, design field 300 is preferably square, even though the shape of 
most lenses is generally round. Defining the field in this shape rather than the usual round 
shape permits the use of a more efficient element. Fig. 4 shows a representative portion 400 
of design field 300 divided by a grid of horizontal, vertical, and diagonal lines 401. (Dashed 
line 400 in Fig. 3 indicates the approximate relative size and location of portion 400 with 
respect to field 300; for clarity of exposition, the vertical lines / and the horizontal lines j are 
here taken to form squares 402 that are 1mm on a side, although normally they are somewhat 
larger. The entire field 300 is divided into elements in a similar manner. Equations 410 
define these lines on the field, assuming a center origin and normalization to unit length. For 
a typical 80mm diameter ophthalmic lens, a typical grid element size (i.e., spacing between 
adjacent horizontal or vertical lines 401) may be in the range of about K4mm. A field 300 
triangulated in this manner may be described as having a uniform triangulation of type II. 

Block 232 introduces a set of basis fimctions for solving the linearized partial 
differential equation. A triangulation of type II permits the basis fimctions to comprise a 
particular type of spline fiinction having compact support and including piecewise quadratic 
polynomials which are continuously differentiable. Employing quadratic polynomials, rather 
than those of higher degree, considerably eases the computational burden over previous 
approaches. 

Fig. 5A graphs a representative B-spline fimction 500 which is piecewise quadratic 
over the octagon 5 10 in Fig. 5B. The numerals at the vertices in Fig. 5 A show the normalized 
coordinates relative to the center [0,0]. This example quadratic B-spline is defined over the 
coordinate range in Fig. 5B as: 
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B{x, y) = B' (x, y) + (x, y) ^B'(x,y)^ B^ (x, y) + B'' (x, y) 

The spline is defined as identically zero outside octagon 510, giving it a compact support. 
The constituent terms of the spline are defined as: 
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This example of a piece wise quadratic spline function is described in P.B. Zwart, 
"Multivariate Splines with Nondegenerate Partitions," SIAMJ, Numer. Anal 10 (1973), 
pp.665-673. (B-spline function B(x,y) has no relation to bilinear function 5(4 v). The 
difference will be apparent fi-om the context.) 
1 0 The generic spline Bix,y) is copied once for every square 402, Fig. 4, over the entire 

field 300 of Fig. 3, and for one layer outside field 300 on all four sides, to form an extended 
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field. If the grid on field 300 is normalized to a unit square [0,l]x[0,l] and has m rows and n 
columns, the spline copies are denoted as: 

for /=0,...,m+l and ;=0,...,«+l. The center of the /j-th spline lies at the point (j- 
5 VjVn]. The spline coordinates may then be scaled to fit the grid selected for a particular 

design. 

The square elements 5 11 of support 5 10 are the same size as square elements 402 of 
field 300. The central square 512 of spline 500 is placed over each square 402 of the 
extended field 300. That is, the outer rows and columns of the support extend one grid 

10 spacing beyond the edge of design field 300 in every direction. Using the splines as basis 

functions realizes the solution to the variational equation as a linear superposition of a number 
of these overlapping splines. That is, a computer solves the equation B{6,v)-L{6) as a set of 
simultaneous equations. First, we set v(x,y)= Hc^jB^jix^y) over all ij\ Substituting this into 
5(<5,v)=L(5) yields a system of linear simultaneous equations with the Cy as unknowns. 

1 5 Solving this linear system gives the values of the Cy, which are then inserted back into 

Hc^jBijix.y) to find the values of v{x,y) for all values of x,y. 

Block 233 imposes the boundary conditions on the linearized equation B{d,v)=L{d), 
We consider three types of boundaries: clamped, partially clamped, and natural (or free) 
boundary conditions; firee and partially clamped conditions may also be grouped together as 

20 less than fiiUy clamped conditions. A clamped boundary specifies that w=0 and (m^,m^)=0 over 

the entire boundary. A partially clamped boundary specifies w=0 and {u^,Uy)=0 for a 
designated part of the boundary. A free boundary allows «, u^, and Uy to have any values on 
the boundary. Our methods permit a unique solution for all three types of boundaries, as long 
as a>(XQ>0 and yS^O. The desired boundary conditions are imposed by selecting only certain 

25 of the simultaneous equations and/or by setting certain variables to desired values. For many 

ophthalmic applications, we have found that natural (firee) boundary conditions produce very 
good results. 
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In this case, using all of the splines overdetermines the solution, so a solution might not exist. 
Dropping the splines in the four comers of design field 300 offers a convenient way to 
achieve a unique solution for natural boundary conditions. This is equivalent to specifying the 
value of v(x,y) at each comer of field 300. 

Should partially clamped conditions be desired for other applications, a designer may, 
for example, set v(jc,l)=0 to clamp the upper edge of design field 300 to the far- vision 
spherical surface. (y=l here indicates the upper edge of a design field normalized to a unit 
square.) In that case, the top row of splines is discarded; that is, those splines whose support 
exist only within the top surface participate in the solution. Fully clamped boundaries may be 
imposed in a similar manner. We believe, however, that a fully clamped boundary produces 
less desirable designs for many ophthalmic purposes. 

Finally, block 234 adds the variations v(x,y) to the base-radius spherical surface u^ix^y) 
to produce a set of heights w=Wo+v above a plane at every point within design field 300» 
(These heights may be easily converted to any other form desired for device fabrication or 
other purposes.) Because the splines exist throughout the design field, and not only at the grid 
points of the triangulation, a computer may calculate and output them for any set of points 
desired by a grinder, plastic casting apparatus, or other tool that actually fabricates a lens. The 
data may be tmncated or discarded outside the edge 100, Figs. 1 and 3, of the desired lens 
diameter. 

Block 240 evaluates the design produced in block 230. An example of such an 
evaluation is to display or print images such as 100 and 1 10, Figs. 1 A and IB for an operator 
to inspect. The spline functions employed in the described embodiment of block 230 are 
quadratic, in order to reduce the amount of computation. However, such splines and their 
superpositions do not have well-defined second derivatives at the grid points. Therefore, it is 
not in general possible to calculate the lens power or astigmatism at some points, because the 
curvature cannot be determined there. Block 240 may then interpolate between the grid points 
with a set of higher-degree tensor-product splines. Fifth-degree tensor-product splines, for 
example, have well-defined second derivatives at all points, and can match the second 
derivatives of adjacent splines at their ends, so that curvature and astigmatism is well-defined 
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at every point in design field 300. The evaluation data are thus calculated from the higher- 
degree splines. Using splines of degree higher than the quadratic splines of block 232 for 
evaluation does not significantly increase overall computation time; the computational effort 
for evaluation is generally much less than for the design. Other evaluation methods, manual 
or automatic, are also possible. 

Blocks 250 output previously calculated data. Block 241 may present the evaluation 
data from block 240 on a display or printer for inspection by the designer. Block 242 outputs 
the design data in a desired format. This format may comprise a data set for direct fabrication 
of a lens according to the design. Outputting the evaluation and design data may involve 
storing or communicating the data, or manipulating it in other ways. For many applications, 
method 200 is may operate on a one-pass basis, where the evaluation data is primarily 
employed for comparison with other designs. However, if desired, method 200 can easily be 
configured to repeat blocks 220-252 for successive iterations of the same design, with or 
without new parameters from the designer. If desired, the computer itself may analyze the 
evaluation data, and specify new design parameters for further iterations. 

Block 260 represents manufacturing a physical lens from the data output by block 252. 
This operation may employ any fabrication technique, such as numerically controlled grinding 
and polishing of glass blanks or molding of cast plastic. 

Fig. 6 illustrates an example of a computerized system 700 for hosting the invention. 
Bus 710 connects a number of components, such as a processor 720, internal storage 730, and 
input/output devices 740. Blocks 740 may include input devices 741 such as keyboard and 
mouse, and output devices 742 such as a display and printer. Storage 750 may accept an 
external medium such as 751 for holding instructions and data for carrying out methods 
according to the invention, for holding data for designs of the optical devices produced by the 
invention, and/or for other data, such as input parameters. One or more interfaces 760 
communicate between components 710-750 and external computing facilities, such as server 
770 and lens fabricating machinery 780. Server 770 may constitute a remote facility for 
holding design programs and data, and/or other items pertaining to system 700. Block 780 
may comprise a numerically controlled (NC) grinding engine of any type, molding apparatus 
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for cast plastic lenses, or any similar apparatus. It may accept design data from system 100, 
by communication over a network 761 or from removable medium 751, to produce PAL 
lenses or other optical devices 782 from blanks 781. In some cases, it may be advantageous to 
combine machinery 780 physically with some or all of the components 710-750, receiving 
data directly rather than communicating over an external network 761. Machinery 780 may 
also supply data back to other components for design verification or other purposes. 

Astigmatism Correction 

The basic design methods described above attempt to reduce all astigmatism in the 
lens to zero. However, many ophthalmic lenses include a prescribed amount of astigmatism 
or cylinder in order to correct a measured amount of astigmatism in a patient's eyes. Single- 
power lenses commonly employ a spherical (or aspherical) surface, and a toric or atoric 
surface to give both power and astigmatism correction in the same surface of the lens. 

In conventional designs, the complex surfaces of progressive lenses preclude 
astigmatism correction in the same surface as power correction. Progressive lens designs, 
including those using the basic methods above, grind or mold the progressive magnifying 
powers into one surface of a lens (usually the front), and fabricate the prescribed astigmatism 
or cylinder correction with toric curves in the other surface, because the power of a thin lens is 
always the difference in powers between the two surfaces. Disadvantages of this technique 
include added time and expense in fabricating both lens surfaces, and optical errors arising 
from the distance between the two surfaces. Further, researchers in recent years have found 
that placing progressive power correction on the back surface of a lens offers advantages such 
as increasing the area of clear vision; placing prescribed astigmatism on the back surface 
denies this feature to patients who require astigmatism correction. 

An extension to the basic methods of the present invention, however, permit a 
prescribed amount of astigmatism correction to be designed and fabricated into the same 
surface as the progressive power correction. The basic methods design progressive powers as 
perturbations from a base surface such as a sphere, or even a plane. However, the equation 
forms chosen for the basic methods do not significantly limit the shape of the base surface. 
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Therefore, the extended methods place a prescribed astigmatism correction into the base 
surface, and design the progressive power correction as a perturbation to that surface. The 
base surface may constitute a surface that also incorporates part or all of the far-view 
correction for the progressive prescription. Although the following examples employ toric 
base surfaces, atoric surfaces are also suitable. 

A torus is a volume of revolution, part of whose surface can be expressed as: 



In this example, the astigmatism correction is horizontal (angle n/2) if R>Ri. Rotating the 
sxirface produces corrections at any other prescribed angle. 

Following the same approach as in the basic methods, we assume that the final surface 
is a base surface with a perturbation, u{x,y)=UQ{x,y)-^v{x,y). In this case, however, the base 
surface is a toric surface that includes at least the astigmatism correction, and may include 
some power correction as well That is, Uq=z on a rotation of z. 

The expressions for ^(m) and K{u) in the previous section could have assumed several 
different forms; the particular forms were selected in order to facilitate the present purpose as 
well. Here again, the error terms and E2 are insignificant as long as x^-^y^^cR^, where c<l. 
The expressions then become: 



w being a dummy variable. We now represent the cylinder — ^i,e., the measure of its 
astigmatism — ^as A(u)-A(uq)-^A^(vX where the first term represents the cylinder of the original 
surface Uq and the second term is the change in cylinder resulting from a perturbation. 




K{u)« K^(u,)+ K^(v)+ K^(u,,v) 



where: 
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In contrast to the basic methods described in the previous section, we desire to keep 
the prescribed amount of cylinder of the original surface, at least in the critical regions 321- 
323, Fig. 3, for far, near, and intermediate vision. Therefore, the goal is to minimize y4„o(v)^ 
instead of A{u)\ Clearly, A{uof ^H^oiuof-K^{uo\ and approximately A^{vy^H^{vf-K^{v), 
We then define the following functional, integrating as usual over the entire design domain 
300, Fig. 3: 

I(u)= l[a{x,y){H^{vf -K^{v))+ P{x,y){H{u)- P{x,y)f]dxdy 

Assume that the surface u{x,y) minimizes the fimctional. As before, let the test function d{x,y) 
be infinitely differentiable over the continuum. Calculating the first variation of /(w) at u and 
using the notations ^^„o(^) K^d^.v) as in the previous section yields: 

\{2{a + p)H^{v)H^{d)-aK^{d,v)]dxdy 

This result yields bilinear and linear terms that are similar to those in the equations for 
B{6,v) and L{d) in the preceding section, the only difference residing in the second line. A 
spherical surface is a special case of a toric surface, and it can indeed be shown that these two 
equations match when is a spherical surface of radius R, 

Method 200, Fig. 2, can now be adapted to design progressive lenses with astigmatism 
correction on the same surface. 

Block 220 adds another input parameter, a conventional specification of a prescribed 
astigmatism power (cylinder) and axis angle. From this, method 200 may easily calculate an 
appropriate toric base surface UQ{x,y), 

Block 231 solves an equation B' {6,v)^L\d) in the form written directly above, rather 
than in the form described in the previous section. This equation is already linearized, 
therefore greatly reducing the computing time and power required for its solution. The 
solution still minimizes the astigmatism toward zero, using the weight function cc, because 
this solution includes only the undesired astigmatism from the progressive powers, and does 
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not carry any part of the prescribed correctional astigmatism. The splines 232 of Fig, 5 or 
others may serve as basis functions as in the previous section. Boundary conditions 233 are 
employed in the same way. 

Operation 234 adds the perturbations from the toric base surface UQ{x,y) defined in 
modified block 220 to produce a set of heights u^Uq^v above a plane at every point within 
design field 300. Again, these heights may be converted to other units and formats if desired. 

The output design data 252 now represents a single surface incorporating both 
progressive power corrections and astigmatism corrections. Block 260 thus need only 
fabricate a single lens surface for both corrections. 

Fig. 6 shows an example ophthalmic progressive addition lens (PAL) 600 that 
incorporates both multiple magnifying powers and prescribed astigmatism in a single surface 
according to the invention. This lens, having an edge 610, typifies designs made with our 
procedure, although other optical devices can be produced in a similar manner. The 
coordinate axes denote millimeters from the center of the lens. Fig. 6A includes contour lines 
620 having constant optical power. Fig. 6B includes contours 630 of constant astigmatism or. 
cylinder. Short bars 640 indicate the direction of maximum principal curvature, that is, the 
direction of the total astigmatism, both for the prescribed correction and for errors induced by 
the progressive powers. 

This example illustrates a lens design of radius 30mm, having an add power of 2,0 
diopters, and a prescribed (i.e., desired) astigmatism correction of 3.0 diopters at an angle of 
7i/4.(northwest in Fig. 6A and B). This design achieves the add power, and the total 
astigmatism is close to the prescribed value in the far- vision, near- vision, and corridor 
regions. The direction of the total astigmatism remains close to the prescribed direction, and 
both the power and the cylinder vary smoothly, without visibly annoying sharp changes. 
There is a clear corridor connecting the far and near regions in the cylinder contour in Fig. 6B. 
The maximum deviation in total astigmatism from the prescribed value is only 75% of the add 
power, that is, 75% of the difference between the magnifying powers of the far- and near- 
vision regions. This is a remarkable result in comparison to previous progressive designs 
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known to us. Typically, a design is considered good when its total astigmatism deviates from 
the prescribed value by less than 100% of the add power over the lens surface. 

System and Products 

Fig. 7 illustrates an example of a system 700 for hosting embodiments of the 
invention. A bus 710 interconnects a digital processor 720, internal storage such as RAM 
and/or ROM 730, and input/output equipment 740, including display and/or printer devices 
741 and keyboard and/or pointing devices 742. I/O equipment 740 may receive input 
parameters for methods according to the invention, and display prompts, evaluations, and 
other data for the methods. Bus 710 also couples external storage devices 750 such as 
magnetic and/or optical drives capable of reading and/or writing removable media such as 
disk 75 1 . One or more interfaces 760 mediate data transfer among bus 7 10 and other 
equipment such as remote server 770 and/or fabrication machinery such as lens grinder 780. 
One or more links 661 may be realized as direct wire or wireless connections, local or wide- 
area networks, or public networks such as the Internet. 

Instructions and data for practicing methods according to the invention may be held 
and/or transferred by various media, such as external storage 651, communication links 661, 
and storage 771 associated with a remote server. 

Numeral 790 illustrates an optical product according to the invention, in this example 
ophthalmic spectacles or eyewear. Machinery 680 accepts lens blanks 791 or other raw 
materials, fabricating therefrom finished progressive lenses 792, possibly with prescribed 
astigmatism correction. Fabricating lenses 792 may include cutting them to a shape different 
from that of blanks 791. Lenses 792 each have a front surface 793 away from the wearer, and 
a back surface 794 toward the wearer. 

If astigmatism correction is included, the basic methods described above may place the 
multiple magnifying powers on one surface such as front surface 793, and the prescribed 
astigmatism on a cylindrical, toric, or atoric back surface 794. Extended methods such as 
those described in the Astigmatism Correction section are capable of grinding, molding, or 
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otherwise fabricating both the progressive magnifying powers and the prescribed astigmatism 
correction directly in a single siuface of blanks or materials 791, in a single operation. 

Some researchers believe that placing the progressive powers on back surface 794 has 
advantages, such as a larger clear- vision area. Methods according the present invention 
permit this where astigmatism correction is also included, by integrating both the magnifying- 
power corrections and the astigmatism correction into the same surface. Placing both on the 
same surface also eliminates undesirable offsets between the two forms of correction when 
they are placed on physically separated surfaces. 

Where astigmatism correction is not included, or where it is included on the same 
surface as the magnifying corrections, the other surface of lenses 792 may have a constant 
power, such as a spherical surface. The two lenses may of course differ from each other in 
magnifying and/or astigmatism corrections, as a particular wearer may require different 
corrections in each eye. 

Spectacles 790 also include any type of frame 795, attached to lenses 792 by 
conventional means such as fasteners or rims (not shown). 

Conclusion 

The foregoing description and the drawings illustrate specific embodiments of the 
invention sufficiently to enable those skilled in the art to practice the invention in those 
embodiments and in other embodiments. Such other embodiments may incorporate structural, 
logical, electrical, process, and other changes. Examples described herein merely typify 
possible variations. Individual components and fixnctions are optional unless explicitly 
required, and operations may be performed in any order, including in parallel Portions and 
features of some embodiments may be included in or substituted for those of others. The 
scope of the invention encompasses the fiiU ambit of the claims set forth below, and all 
available equivalents. 

We claim as our invention: 



-20- 



